Josephson current in strongly correlated double quantum dots 
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We study the transport properties of a serial double quantum dot (DQD) coupled to two superconducting 
leads, focusing on the Josephson current through the DQD and the associated 0-7T transitions which result 
from the subtle interplay between the superconductivity, the Kondo physics, and the inter-dot superexchange 
interaction. We examine the competition between the superconductivity and the Kondo physics by tuning the 
relative strength A/TV of the superconducting gap A and the Kondo temperature TV, for different strengths of 
the superexchange coupling determined by the interdot tunneling t relative to the dot level broadening Y. We 
find strong renormalization of t, a significant role of the superexchange coupling J, and a rich phase diagram of 
the and tt -junction regimes. In particular, when both the superconductivity and the exchange interaction are 
in close competion with the Kondo physics (A ~ J ~ TV), there appears an island of 7r'-phase at large values 
of the superconducting phase difference. 

PACS numbers: 73.23.-b,72.15.Qm,74.45.+c 



Introduction. — In a metal containing a dilute concentra- 
tion of magnetic impurities the competititon between Kondo 
physics, which favours screening of the localized spins by 
the itinerant conduction-band electrons fl|], and antiferromag- 
netic (AF) exchange interactions between impurities leads to 
a quantum phase transition |2]. Even more interesting prop- 
erties emerge when the metal turns superconducting below 
the critical temperature. For s-wave superconductors, Cooper 
pairs formed by itinerant electrons |3J are yet another possi- 
ble singlet state which competes with the above. The intrigu- 
ing interplay of these phenomena, which might actually co- 
exist in complex materials such as heavy-fermion supercon- 
ductors, governs the low temperature physics of these sys- 
tems. Nanoscale systems allow to tune the ratio between the 
relevant parameters (the Kondo temperature TV, the AF ex- 
change interaction J, and the superconducting gap A, respec- 
tively) and, therefore, enable thorough investigations of such 
competition in a controlled setting. In the simplest case of 
single quantum dots attached to superconducting reservoirs, 
where only Kondo physics and superconductivity are rele- 
vant, a sign change of the Josephson current, from positive 
0-junction to negative tt -junction behavior, signals a quantum 
phase transition between a singlet and a doublet ground state 
as TV// A decreases J3-01- This to tt -junction transition has 
been experimentally realized, confirming some of these phys- 
ical aspects (7)]. A double quantum dot (DQD) coupled to 
normal metals constitutes a physical realization of the two- 
impurity Kondo model 181 41011 . as demonstrated experimen- 
tally III lUl 211 . When the reservoirs become superconducting, 
this system is a minimal artificial realization of the described 
competition among three different spin-singlet ground states. 
In this Letter we focus on a detailed analysis of the Josephson 
current which, as a ground state property, shows signatures of 
this subtle competition. Our results, obtained by a highly ef- 
ficient numerical renormalization group (NRG) scheme, able 



to deal with superconducting correlations, are summarized in 
Figs. |2] and [3] The interplay between superconductivity and 
the Kondo physics is studied by tuning the relative strength 
A/TV. The role of the superexchange coupling is tuned by 
the interdot tunneling t relative to the dot level broadening 
r. Importantly, we find strong renormalizations of t and a 
significant role of the superexchange coupling J compared 
with the previous works H, @, 13] based on the slave-boson 
mean-field theory (SBMFT). Moreover, we find a rich phase 
diagram of the 0-ir transition. In particular, when both the 
superconductivity and the superexchange are in close com- 
petion with the Kondo physics (A ~ J ~ Tk), there appears 
an unexpected island of 7r'-phase at large values of the su- 
perconducting phase difference (f>. We provide clear interpre- 
tations by examining the spin-state- resolved Andreev bound 
states inside the superconducting gap. 

Model. — The system that we consider is a DQD modelled 
as a two-impurity Anderson model connected to two super- 
conducting leads (Fig.[TJi described by standard BCS Hamil- 
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FIG. 1: (color online) Schematics of the DQD system coupled to 
superconducting leads. In the deep Kondo limit, this system is an 
artificial realization of the two-impurity Kondo problem in the pres- 
ence of superconducting correlations. 
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tonians 13]: % = Hd + %l + Ut, where 

(2) 
(3) 
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Here describes an electron with energy e^, momentum 
k, and spin p on the lead £ = 1,2, and an electron in 

the dot i = 1,2; n £k =J2fj. c lk^ c ^ and 71 » = S M 
e is the single-particle energy on each dot that is tuned by 
gate voltages and U is the on-site Coulomb interaction. The 
electrons can tunnel between the two dots with the tunneling 
amplitude t. is the superconducting gap, and <j>i the phase 
of the superconducting order parameter. 

The two leads are assumed to be identical except for the 
superconducting phases. Assuming that the leads have a flat 
band with the density of states p = 1/2D, where 2D is the 
bandwidth, the hybridization between the dots and the leads is 
well characterized by a single parameter T = irpV 2 . Since we 
are interested in the Kondo correlations, we concentrate on the 
Kondo regime with localized level — e ^> T and large charging 
energy U > 2|e|. For the representative results shown below, 
we choose T = 0.014D or T = 0.02D, fix e = -0.2D, and 
take the large U — oo limit. The parameter space is never- 
theless still large. We examine the results by varying A/Tk, 
t/T, and (f> = (f> L - (f> R . 

We solve the Hamiltonian using the numerical renormaliza- 
tion group (NRG) method 114]. Due to the relatively low sym- 
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conservation), the NRG iteration is numerically very demand- 
ing and special attention is necessary to obtain reliable results. 
Using a new discretization scheme the numerical artifacts due 
to a large discretization parameter A are almost completely 
canceled out; the technical details are given elsewhere. 

We first review two crucial effects from a previous work on 
DQD with normal leads fioll which keep playing significant 
roles in the present superconducting case. Firstly, the inter- 
dot tunneling t is significantly renormalized compared with 
the predictions based on the SBMFT 12]- It is impor- 



tant to treat this effect properly to fully account for the trans- 
port. Secondly, the inter-dot antiferromagnetic superexchange 
J is finite even for U = oo. lfioll This interaction is mediated 
by the virtual tunneling of electrons to the conduction leads. 
It is thus important to take into account the interplay of the 
superexchange coupling with the superconductivity and the 
Kondo effect. 

Strong coupling limit (Tk 3> A). — Figure [2] shows nor- 
mal state conductance (a,b) and the critical Josephson current 
(c,d). Here we focus on the strong coupling limit, A/Tk = 
0.1 (black circles). In this limit, the critical Josephson cur- 
rent / shows similar features as the normal-state conductance. 
It peaks at equal values of t/T and, remarkably, when the 
conductance in the normal state is unitary, the Josephson cur- 
rent reaches the quantum limit If = eA/h. This is expected 
since the Kondo effect dominates over the superconductivity, 
therefore the transport is determined by the competition be- 
tween the Kondo physics and the interdot superexchange (for 
t/T < 5) or interdot molecular orbitals (t/T > 5). As we 
analyzed in detail in the previous work for the normal lead 
case, the peaks in G and I c at t/T « 0.4 (for T = 0.014D) 
result from the crossover from the "Kondo singlet" to the "su- 
perexchange singlet": For t/T < 0.4, J < Tk, wheares for 

0. 4 < t/T < 5, J > T K . We stress that the crossover is 
shifted significantly to smaller t/T = 0.4 compared with the 
estimation from SBMFT llii l9l [l3ll due to the previously men- 
tioned strong renormalization of t. As t/T increases beyond 
5, the DQD starts to form molecular orbitals and effectively 
behaves as a single QD. In this regime, the Josephson current 
in single dots has n shift in the Coulomb blockade regime 

a- 

01, which is directly related to the 7r-phase in Fig. [3](a,b) for 
large t/T. 

For comparison, we have also calculated the critical Joseph- 
son current as / = max^ I(<f), where 1U5I1 
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FIG. 2: (color online) Normal state conductance (a,b) and the critical 
Josephson current in the superconducting state (c,d) as a function of 
t/T for T = 0.014D (a,c) and T = 0.02D (b,d). The results are 
expressed in units of the conductance quantum Go = 2e 2 /h and 
the supercurrent quantum 7^ = eA/h. In the superconducting case 
(c,d), the different curves are for A/Tk =0.1 (black circle), 3.6 
(red square) , and 10 (blue diamond). The dashed line is from an 
effective non-interacting theory, Eq. I0. 



where 7* = eA/h is the critical current of a transparent 
single-mode junction, and g = G/Gq is the (dimensionless) 
normal-state conductance obtained from a NRG calculation. 
These relations are applicable only under the assumption that 
the QD state is weakly affected by the superconductivity. As 
can be seen in Fig. [2] for t/T < 5, the result from the effec- 
tive theory (dashed lines) and the fulll NRG calculation show 
a qualitative agreement. For larger t (t/T > 5), we enter the 
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FIG. 3: Phase boundaries between the 0- and 7r-states of the Joseph- 
son current for T = 0.014D (left panels) and T = 0.02D (right pan- 
els). When both the inter-dot superexchange coupling and the super- 
conductivity are in close competition with the Kondo effect between 
the superconducting leads and adjacent dots (t/T ~ 0.2, A ~ Tk), 
there appears an island of 7r'-phase at larger phase difference. 



single-dot regime and the Kondo effect is suppressed. This 
leads to the deviation of the two results. 

Weak coupling limit (Tk <C A). — The results, shown in 
Figs. |2] and [3] for A/Tk = 10 (gray diamonds, blue online), 
indicate that the superconducting correlations in the leads sup- 
press the Kondo effect and the Josephson current remains 
small until the system enters the single-dot regime. It is re- 
markable that the S-DQD-S system behaves as a 0-junction in 
the weak coupling limit in contrast to the S-QD-S case where 
the 7T -junction appears in the same regime. In single quan- 
tum dots, the appearance of 7r-junction is due to the rever- 
sal of the order of the electrons forming Cooper pairs after 
tunneling. J3l However, in series coupled DQD, the order is 
preserved so that no additional phase factor arises from the 
tunneling and thus 0-junction is formed even in the Coulomb 
blockade regime. Specifically, the perturbation theory applied 
in the weak coupling limit for U —> oo gives 
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with _Ek = -\/ e k + O ur NRG calculations confirm that 
the ground state is a spin singlet as long as the 0-junction is 
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FIG. 4: (color on-line) Left panels: Josephson current vs 
phase difference near the 7r'-phase (A/Tk = 3.6 and t/T = 
0.1, 0.2, 0.3, 0.4, from top to bottom). Right panels: Corresponding 
spin-state-resolved Andreev bound states inside the superconducting 
gap. Spin singlet states are depicted by (red) circles and doublets by 
(black) squares. The changes in the spin states are closely related to 
the island of 7r' phase in Fig. [3] (c) [similar spin-state-resolved An- 
dreev states for corresponding parameters will also explain the island 
in Fig.[3](d)]. 



formed. Hence, in contrast to the single quantum dot system, 
in a large part of the parameter space there exists no phase 
transition as we move from the weak to the strong coupling 
limits by varying A/Tk- we always have a spin singlet state 
with 0-junction behavior. 

Subsequent transition into ir -junction for very large t/T 
in Fig. [3] (e,f) is ascribed to the competition between effec- 
tive spin- 1/2 Kondo correlation and superconductivity as in 
the strong coupling limit. Since the superconducting gap 
(A/T K 3> 1) is larger than in the strong coupling limit 
(A/T K <C 1), however, the transition can take place at 
smaller values of t/T for which the effective Kondo tempera- 
ture Tk is higher. 

Moreover, the critical current is relatively large even though 
the system is in the weak coupling limit, while in the sin- 
gle quantum dot the critical current in this limit is very small 
(I / 1 & C < 0.1) |@]. A very likely explanation is that the one- 
electron spin- 1/2 Kondo state is formed at smaller values of 
t/T and that strong superconductivity is responsible for it. 
The (one-electron) Kondo assisted tunneling then makes the 
junction more transparent and enhances the critical current. 
Hence, the physical origin of the peak in the critical current is 
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different in the weak and strong coupling limits. 

Intermediate coupling (Tk ~ A). — A highly non-trivial 
behavior occurs for Tk ~ A. In this regime, the superconduc- 
tivity, the superexchange, and the Kondo physics can all be in 
close competion. This subtle interplay keeps the Josephson 
critical current finite, somewhere between the weak coupling 
and strong coupling limit, Fig. [2] (c and d), except for very 
large t/T, where single-dot physics again governs the trans- 
port. 

More interestingly, the phase diagrams in Fig.[3](c,d) reveal 
the reentrance behavior from the 0-j unction, to tt' -junction, 
and back to 0-junction, and eventually to the tt -junction for 
larger superconducting phase difference <\>. In order to under- 
stand this behavior, we closely examine the subgap Andreev 
bound states, which are Bogoliubov quasi-particle excitations 
from the ground stated and whose derivatives with respect to 
<t> give the Josephson current lfl6ll . In Fig.|4]we plot the energy 
levels of the Andreev states as a function of <p, in the param- 
eter regime corresponding to the island in Fig. [3(c). Let us 
focus on, say, </> = 0.97T. For t/T < 0.1, the singlet Kondo 
clouds are formed between the superconductors and the adja- 
cent dots, thus the ground state is likewise a spin singlet, while 
the excitations correspond to doublet states, see Fig.|4](b). As 
the Josephson current is given approximately by the phase- 
difference-derivative of the Andreev levels, this corresponds 
to the 0-junction behavior; Fig.@](a). For 0.4 < t/T < 5, the 
local inter-dot singlet state is induced on the DQD due to the 
antiferomagnetic superexchange interaction, thus the ground 
state is again a spin singlet; Fig.|4](h). As before, this results 
in the 0-junction behavior; Fig. |4](g). 

In the above two cases, both Kondo effect and superex- 
change barely win over the superconductivity, for all values 
of <fi. However, when 0.1 < t/T < 0.4, the Kondo effect is 
suppressed by the large phase difference. Q This is indicated 
in the fact that the ground state is now a doublet, while the ex- 
cited state is a singlet, as shown in Fig. H](d,f). Accordingly, 
the transport properties are different and, in particular, the tt- 
junction behavior is observed; see Fig. |4](c,e). This regime 
is denoted as n' in the phase diagram in Fig. [3] While the tt 
phase for large t corresponds to the single occupancy of the 
dots, the 7r' phase occurs for the double occupancy. For large 
values of T (T = 0.02ZT), the tt' island becomes bigger and it 
can merge with the tt regime. 

Conclusion. — We have studied the Josephson current 
through a serial double quantum dot coupled to two supercon- 
ducting leads. We have observed a strong renormalizations 
of t, the significant role of the superexchange coupling J, 
and a rich phase diagram featuring different 0-tt transitions. 
In particular, when both the superconductivity and the su- 
perexchange were in close competion with the Kondo physics 
(A ~ J ~ Tk), there appeared an island of 7r'-phase at larger 
values of the superconducting phase difference. This finding 
motivates further studies of this regime, which may shed new 
light on the physics of heavy-fermion superconductors. 
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